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/**\ , The dosimetry for radiocolloid therapy of cystic craniopharyngiomas is investigated. Analytical 

calculations based on the Loevinger and the Berger formulae for electrons and photons, respectively, 
are compared with Monte Carlo simulations. The role of the material of which the colloid introduced 
inside the craniopharyngioma is made of as well as that forming the cyst wall is analyzed. It 
is found that the analytical approaches provide a very good description of the simulated data in 
the conditions where they can be applied (i.e., in the case of a uniform and infinite homogeneous 
medium). However, the consideration of the different materials and interfaces produces a strong 
reduction of the dose delivered to the cyst wall in relation to that predicted by the Loevinger and 
'■^ ■ the Berger formulae. 
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I. INTRODUCTION 

C/2 ■ 

O ■ 

' 55 ' Craniopharyngiomas are tumors showing an incidence of 3% for all tumors in adults and 6 to 10% of tumors in 
^*-y children [l|. Though histologically benign, they are effectively malignant because they usually appear in a situation 
' which can affect to important organs such as hypothalamus, optic nerves or chiasms. 

Treatment based on surgery followed (or not) by radiation therapy needs a total resection to guarantee a good 
chance of cure. However, a complete excision is rarely possible due to the vicinity of the organs mentioned above y|. 
7-H ! Alternatively, the introduction of radioactive colloids into the cyst has been considered since the early 50 's 0.O|. 
!> ■ Cyst wall constitutes the target volume and, as a consequence, /3-emitters, as e.g. 32 P and 90 Y are the most adequate. 
Also /37-emitters, as 186 Re and 198 Au, have been used, since 7-radiation allows to test the colloid distribution by 
*^ obtaining a gamma camera image. 

Cystic craniopharyngiomas are roughly spherical and the wall of the cyst presents a thickness between 1 and 3 mm. 
Once the radiocolloid is introduced inside the craniopharyngioma, part of it is attached to the inner surface of the 
£f*l • cyst, while the rest appears to be distributed into the inner volume. 

Dosimetry calculations for (3 emitters have been carried out either on the base of the Loevinger formulation Q (see 
e.g. @) or by using the (3 point kernels generated by Berger (see e.g.||). In the case of /^-emitters, the dose 
rates for photons have been evaluated following the approach of Berger |9| in terms of the build-up factors for the 
"^5 ■ appropriate energies (see e.g. @). 

In order to perform the dosimetry, because of the characteristic distribution of the radiocolloid, two extreme 
*£h situations have been considered in practice [y, Q • The first one is the spherical shell source, where the emitter is 
Oh 1 supposed to be uniformly distributed in the inner surface of the cyst. The second one is the spherical volume source 
, (SVS), which assumes the radionuclide uniformly distributed in the inner volume of the craniopharyngioma. In all 
r *~J ■ cases, a uniform and infinite homogeneous medium is considered and the presence of different material media and 
rS ' interfaces is not taken into account. 

Monte Carlo (MC) simulation offers the possibility to perform a more realistic dosimetry. In this paper we have 
studied how the consideration of the actual media involved in the craniopharyngiomas modifies the results of the 
standard calculations. 



II. MATERIAL AND METHODS 



Here we have assumed that the cystic craniopharyngioma is described by means of a spherical layer of thickness 
S and inner radius R (see Fig. 1). We adopted SVS approach and then we dealt with extended uniform spherical 
sources of given radionuclides. We were interested in calculating the dose delivered to the cyst wall as well as to 
points external to the craniopharyngioma and related to the critical surrounding tissue. The corresponding dose rate 



FIG. 1: Model scheme used to describe the craniopharyngioma. 



is accounted for by evaluating the integral 

D(x) = [ dr C(r) J(x — r)=A [ drJ(|x-r|), (1) 

where V represents the source volume, x and r are, respectively, the positions of the target point and of the volume 
element of the distribution with respect to its center, C(r) = A gives the radionuclide activity concentration and 
J(|x — r|) represents the absorbed dose rate in x due to a point source in r. 



A. sources 



In the case of f3 sources, dosimetry calculations can be performed on the basis of the Loevinger formula |5|, which 
gives the absorbed dose distribution around a point source of a /3-emitting radionuclide in a homogeneous infinite 
medium. In this approach the absorbed dose rate at a distance x from the point source is given by 



Jp{x) = Ti(x) + T 2 (x) , 
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In previous equations, z = (pv) -1 and t — cz are characteristic distances, p the density of the homogeneous medium, 
v an apparent absorption coefficient and c a dimensionless parameter. The normalization constant B is evaluated by 
imposing that all the emitted energy is absorbed in a very large sphere, and is given by 
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and Ef3 the average energy per disintegration. Finally, 

0(r) - 



r < 0; 
r > 0. 



The parameters v and c are characteristic of each /3-emitter. For the first one, Loevinger et al. gave the following 
numerical parameterization: 
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where £'™ ax is the maximum energy of the electrons emitted by the source, in MeV, and v is in cm 2 g _1 . Ep is the 
average energy for a hypothetical (3 disintegration of allowed type with the same maximum energy. The c parameter 
was parameterized as 5] : 



2 , 0.17 MeV < £™ ax < 0.5 MeV : 
<■ = < 1.5, 0.5 MeV < Ef ax < 1.5 MeV; 



1, 



1.5 MeV < Eff** < 3 MeV. 



Adopting this approach for the point absorbed dose rate, the integral in Eq. (JIJ can be evaluated analytically (see 
Appendix) and using Eq. J2<| we can write 
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B. Photon sources 



When /37-emitters are used, the photon dosimetry is usually based on the formula of Berger In this approach, 
a photon point isotropic source is assumed to deposit a dose rate at a distance x from the source which is given by 



J-y(x) = T~ ~ E t ~7 exp(-^x) B cn (fix) . (6) 

47T p X z 



Here fi en and /i are the linear photon energy-absorption and the linear photon attenuation coefficients, respectively, 
at the energy of the emitted photon, _E 7 , and B Qn is the energy-absorption buildup factor, which takes into account 
the contribution of the scattered photons. Following Ref. 0, the buildup factors are expanded as 
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B cn (fix) =J2 b " (rf . (7) 

n=0 

Here bo — 1 and the remaining coefficients were calculated in for some energies ranging from 10 keV to 3 MeV. 

As in the case of the Loevinger formula, the integral in Eq. (JIJ can be calculated analytically for the dose rate in 
Eq. ©, in the case of the uniform spherical radionuclide distribution we are assuming here (see Appendix). From 
Eq. ifTTTT) we have 

D j{x) = -A^-E. < Y J r^K^x), (8) 
p x j=0 

where the different terms are defined in Eqs. I|A12|) - (|A17|) . 



C. Simulation procedure 

MC simulations were performed by using the code PENELOPE [Hj. PENELOPE is a general purpose MC code 
which allows to simulate the coupled electron-photon transport. Analog simulation is performed for photons. For 
electrons, the simulation is carried out in a mixed scheme where collisions are characterized as hard and soft. After 
fixing a value for a critical angle, the collisions with scattering angles larger than the critical value are called hard 
collisions and are simulated individually. The collisions with scattering angle smaller than the critical value are called 
soft collisions and are described by means of a multiple scattering theory. The electron tracking is controlled by four 
parameters. Two of them, called C\ and C2, refer to elastic collisions. The first one, C\, gives the average angular 
deflection due to an elastic hard collision and to the soft collisions previous to it. The second parameter, C2, represents 
the maximum value permitted for the average fractional energy loss in a step. The other two parameters, called W cc 
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and W cr are energy cutoffs to distinguish hard and soft events. Thus, the inelastic electron collisions with energy loss 
W < W cc and the emission of Bremsstrahlung photons with energy W < W cx are considered in the simulation as soft 
interactions. 

For arbitrary materials, PENELOPE can be applied for energies up to 1 GeV and down to few hundred eV, in the 
case of electrons, and 1 keV, for photons. Besides, PENELOPE permits a good description of the particle transport 
at the interfaces and presents a more accurate description of the electron transport at low energies in comparison to 
other general purpose MC codes. These characteristics make PENELOPE to be an useful tool for medical physics 
applications as previous works have pointed out (see e.g. Refs. [l2|-|l4jL Details about the physical processes 
considered can be found in Ref. [TT| . 

In our simulations, the cyst was supposed to be surrounded by an infinitely extended water volume. Both, the 
spherical layer simulating the craniopharyngioma as well as the colloid introduced into the cyst were assumed to be 
made of different materials. Thus, we labelled the dose rates obtained with the simulations as D™, where the subscript 
u> refers to (3 or 7 according to the primary particle emitted and the superscripts v and s indicate, respectively, the 
materials of which the inner volume and the spherical shell are made of. Here we considered four materials: soft tissue 
(t), compact bone (b), gel (g) and water (w). The respective compositions are given in Table [I] 



Composition [atoms/molecule] 




Soft tissue (ICRU) Compact bone (ICRU) 


Gel 


Density [g cm ] 


1.0 


1.85 


1.2914 


H 


0.10117 


0.52790 


0.54697 


C 


0.11100 


0.19247 


0.23524 


N 


0.02600 


0.01603 


0.05393 





0.76183 


0.21311 


0.16156 


Mg 




0.00068 




P 




0.01879 




S 




0.00052 


0.00230 


Ca 




0.03050 





TABLE I: Composition of the different materials used in the MC simulations performed in this paper. 

The point sources representing the radionuclide were assumed to be uniformly distributed inside the craniopharyn- 
gioma. The direction of the emitted particles was supposed to be isotropic around the initial position. For electrons, 
the initial energy was sampled from the corresponding Fermi distributions, taken from Ref. |l6j . For photons, the 
initial energy was sampled according to the relative intensity of the different energies of the corresponding spectra. 

Electrons and photons were simulated for energies above 100 eV and 1 keV, respectively. Below these energies, the 
particles were considered to be locally absorbed. The simulation parameters were fixed to the values: W cc = 5 keV, 
W cr = 1 keV, Ci=C 2 = 0.05. 

Due to the spherical symmetry of the adopted SVS geometry, the dose distributions were supposed to be functions 
of the distance to the center of the distribution, x. The full simulation volume was subdivided into spherical shells with 
radial thickness of 0.5 mm. The energy deposited in each voxel was scored to obtain the corresponding histograms. 

The statistical uncertainties were calculated by scoring both the energy deposited in a voxel and its square for each 
history. The average energy deposited in the k-ih bin (per incident particle) is 

1 N 

t=i 

where N is the number of simulated histories and e^k is the energy deposited by all the particles of the z-th history 
(that is, including the primary particle and all the secondaries it generates). The statistical uncertainty is given by 
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In our calculations, 5 • 10 6 histories were simulated in each run. 
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FIG. 2: Dose rate in water for a spherical source of 32 P with R = 1.75 cm, S = 0.1 cm and an activity concentration of 
1 MBq cm -3 . The distance of the target point to the center of the distribution is labelled x. Dots are for the MC calculation 
(Error bars correspond to la and are smaller than the size of the symbol used). Solid curve gives the result obtained from 
the Loevinger formula while dashed-dotted curve correspond to the Vynckier-Wambersie formula. The lower panel shows the 
same results in semilogarithmic scale. Besides the MC results obtained switching it off Bremsstrahlung are plotted with black 
squares. 



III. RESULTS 



A. Pure f3 radionuclides 



We started by considering a pure j3 emitter as radionuclide in the colloid: the 32 P nucleus. The /3-spectrum of this 
isotope has a maximum energy E™ ax — 1.71 MeV with a mean energy Ep = 0.695 MeV. The decay transition to 32 S 
is of allowed type and the parameters entering in the Loevinger formula are c = 1 and v = 9.18 cm 2 g _1 . Initially we 
analyzed cysts with R — 1.75 cm and 5=1 mm, while the activity concentration of the radiocolloid was supposed to 
be 1 MBq cm~ 3 . 

First, we have assumed that both the inner volume and the layer were made of water. In this situation (uniform and 
infinite homogeneous medium) , the results of the simulation can be compared directly with our analytical calculations 
based on the Loevinger approach (Eq. JSJ) or other analytical approaches such as that of Vynckier and Wambersie 
|17|. In this last method the f3 point kernels generated were more accurate than those in the Loevinger approach. 
The results are shown in the upper panel of Fig. 2, where the full and dashed curves show the results obtained with 
the Loevinger and Vynckier and Wambersie approaches, respectively, and the open squares are the results of our 
MC simulation. In the lower panel, the same results are plotted in semilogarithmic scale in order to emphasize the 
differences at large distances. In Table [H] we show the values of the relative differences calculated as: 



'bj»(xj 


A 




B 






B 



(9) 



where the subscripts A and B stand for the Loevinger (L) and Vynckier and Wambersie (VW) approaches or the 
MC simulation. As we can see, Loevinger and MC results agree rather well up to x ~ 2 cm. Above this point, the 
disagreement is apparent, but the value of the dose rate is, at least, one order of magnitude smaller than the dose 
rate in the cyst wall. Vynckier and Wambersie model provides a slightly better result up to x ~ 2.3 cm. At larger 
distances, the Vynckier and Wambersie approach goes to zero, but there the dose rate is, at least, three orders of 
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magnitude below that in the wall. Similar results are obtained if the (3 point kernels used are those provided by Berger 



A,L A ' B) W [%] 

x [cm] (L,MC) (VW.MC) (VW.L) 



1.775 
1.825 
1.875 
1.925 
1.975 
2.025 
2.125 
2.225 
2.325 
2.425 
2.525 
3.025 
3.525 
4.025 
4.525 
5.025 



-1.1 
2.6 
5.5 
10.0 
17.2 
30.7 
92.4 
258.3 
754.5 
921.5 
386.8 
-93.6 
-99.9 
-100.0 
-100.0 
-100.0 



-2.8 
-2.0 
-2.6 
-3.8 
-5.3 
-5.7 
-5.1 
-17.2 
-68.3 
-100.0 
-100.0 
-100.0 
-100.0 
-100.0 
-100.0 
-100.0 



-1.7 
-4.4 
-7.7 
-12.5 
-19.2 
-27.8 
-50.6 
-76.9 
-96.3 
-100.0 
-100.0 
-100.0 
-100.0 
-100.0 
-100.0 
-100.0 



TABLE II: Relative differences (see Eq. in % between the results obtained for 32 P, in the case of unique medium, with the 
Loevinger (L) and Vynckier and Wambersie (VW) approaches and MC simulation. 



To finish the discussion for a unique medium, we analyzed the role of the Bremsstrahlung radiation. PENELOPE 
includes this mechanism and, in order to estimate its importance we have performed a new simulation by switching 
it off. The results obtained are plotted with black squares in the lower panel of Fig. 2. As we can see, the effect 
of the Bremsstrahlung shows up above ~ 2.4 cm. However, at that distance, the dose deposited is four orders of 
magnitude smaller than the dose delivered to the cyst wall. On the other hand, we can see that the Vynckier and 
Wambersie model does not include Bremsstrahlung at all, while the Loevinger approach takes it into account in a 
somehow average way. Despite this, the conclusion is that these analytical approaches permit a reasonable description 
of the dose rate delivered by this kind of extended sources within a unique medium. 

The next step was to investigate the role played by the interfaces present in the craniopharyngioma. The idea 
was to elucidate, by using MC simulation, if the various interfaces and materials modify the dosimetry provided by 
the Loevinger approach. To measure the effect of the presence of the different materials, we calculated the relative 
differences 

A vV-vs (cyst) = ^ V (cy S t)-^(cy S t) ) 

L>^(cyst) 

where 

£>: s (cyst) = (11) 

are the integral dose rates delivered inside the cyst wall. Here k runs over all the bins contained inside the cyst wall, 
p s is the density of the material forming this wall and v s is its volume. 

We have first studied the effect of the material forming the spherical shell which represents the craniopharyngioma. 
We have considered two extreme cases: soft tissue and compact bone. The differences of the dose rates calculated 
with MC in both situations with respect to D™ w (x) are shown in panels (a) and (b) of Fig. 3. The relative differences 
in the cyst wall are given in Table Ifffl As we can see, the presence of the soft tissue in the spherical shell leaves 
practically unmodified the dosimetry. In fact, the total dose rate delivered to this shell increases by only 0.9% in this 
case. On the contrary, when the bone is considered the dose received by the wall reduces 28.0%. Furthermore, the 
dose rate delivered to the external medium nearby the craniopharyngioma diminishes appreciably in the case of the 
bone, while no modifications are observed for soft tissue. Obviously, these modifications are not considered in the 
analytical expressions based on the Loevinger approach. 
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FIG. 3: Effect of the materials considered for the spherical shell which represents the cystic wall in our model and for the inner 
volume. Panels (a) and (b) show the effect of considering the wall of the craniopharyngioma to be done of soft tissue and bone, 
respectively. Panels (c) and (d) show the effect due to the substitution of the water by the gel in the inner volume. Finally, 
panels (e) and (f) show the full effect due to the presence of these materials with respect to the reference dose rate D^f(x). 
The shadow region shows the cystic wall which has been considered to have a thickness 5 = 1 mm. Error bars correspond to 
la and are smaller than the size of the symbol used. 



In the treatment procedure, the radionuclide is introduced in the craniopharyngioma inner volume, embedded in a 
certain type of colloid. In our simulations we have considered the gel with the composition given in Table ||| In panels 
(c) and (d) of Fig. 3 we show the differences in the dose rates obtained when the water is substituted by the gel in 
the inner volume. As we can see, the gel produces a noticeable reduction of the dose delivered to the cyst wall. This 
reduction reaches 20.7% in the case of soft tissue and 20.6% in the case of the bone. Besides, also the dose delivered 
to the surrounding outer volume reduces. 

The net effect due to the presence of both gel in the inner volume and soft tissue or bone in the craniopharyngioma, 
in relation to the case in which only water is considered, is shown in panels (e) and (f) of Fig. 3. The dose rate 
delivered to the cyst wall reduces in 20.0%, in the case of the soft tissue, and 42.9% in the case of bone. Even if we can 
expect that the craniopharyngioma is composed of a material with an intermediate density, the dose rate delivered 
to both, the cyst wall and the outer volume nearby it, would be strongly reduced in relation to what it is obtained 
using the Loevinger approach. 











s '- vs (cyst) [%] 




R [cm] 


S [mm] 


wt-ww 


wb-ww 


gt-wt gb-wb gt-ww 


gb-ww 


1.75 


1.0 


0.87 


-28.02 


-20.67 -20.61 -19.99 


-42.86 


1.75 


2.0 


1.18 


-39.65 


-21.15 -21.24 -20.22 


-52.47 


1.75 


3.0 


0.98 


-44.77 


-20.81 -20.99 -20.03 


-56.37 


1.00 


1.0 


0.46 


-27.39 


-20.62 -20.58 -20.25 


-42.34 


2.50 


1.0 


0.50 


-28.44 


-20.50 -21.02 -20.11 


-43.48 



TABLE III: Relative differences in %, as given by Eq. JTOJ, for 32 P in the different cases analyzed in this paper (see text). 
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FIG. 4: Effect of the Bremsstrahlung radiation for the case in which gel and bone are present in the inner volume and the wall 
of the craniopharyngioma, respectively. Black squares give the dose rate obtained when Bremsstrahlung mechanism is switched 
off. Open squares correspond to the full calculation. 



The presence of materials with densities larger than that of water, could make the Bremsstrahlung effect more 
important. As before, we have carried out a new simulation (for the case in which gel and bone are present in 
the inner volume and in the wall of the craniopharyngioma, respectively) , but we switched off the Bremsstrahlung 
mechanism. In Fig. 4 the results obtained (black squares) are compared to those of the full calculation (open squares) . 
These results are rather similar to those found for the case of a unique medium (see lower panel in Fig. 2), except 
for two details. First, the Bremsstrahlung effect begins to be observed in the radiation queue at a slightly shorter 
distance, ~ 2.2 cm. Second, the inclusion of Bremsstrahlung effects produces an increase of the dose delivered to the 
cyst wall which reaches the non negligible value of ~ 45%. If we compare the simulation without Bremsstrahlung with 
the analytical approaches which, as mentioned previously, do not treat this effect in an adequate way, the reduction 
observed for the dose delivered to the cyst wall would be even larger. 

At this point, it is necessary to clarify how the previous conclusions are related to the choice of the geometry. For 
this reason we have performed a set of calculations similar to the previous ones but assuming different thicknesses for 
the tumor wall and different sizes for the inner volume. 

First we have analyzed the effect of the wall thickness and we have considered 6 = 2 and 3 mm. The relative 
differences A^ s _vs (cyst) obtained are shown in Table 11111 for different materials and configurations. As we can see, 
the effect of the presence of the soft tissue is rather small, increasing the dose delivered to the cyst by ~ 1% (see third 
column). On the contrary, the fact that the wall is made of bone produces a considerable reduction of the dose rate 
in it. Besides, this reduction grows with the thickness of the cyst wall (see fourth column). 

The effect of the presence of the gel in the inner volume seems to be almost independent of the thickness of the 
wall producing a reduction of ~ 20% in all cases (see fifth and sixth columns). 

Finally, the combined effect of soft tissue in the wall and of the gel in the inner volume (seventh column) is practically 
equal to that of the gel alone, while in the case of bone (eight column), the reduction increases with the wall thickness, 
following the behavior observed when no gel was considered. It is worth to point out that the dose rate delivered to 
the cyst can be diminished by more than 50%. 

To study the influence of the inner volume of the craniopharyngioma we have considered two different values of its 
radius: R = 1 and 2.5 cm. In both calculations we have maintained the activity concentration of 1 MBq cm -3 . The 
results obtained in the simulations are shown in the two last rows of Table ITTT1 As we can see, the values of the relative 
differences are very similar to those found for the former calculation (R = 1.75 cm, 6=1 mm). This indicates that 
the important parameter is the thickness of the cyst wall. The reduction of the dose delivered to it becomes bigger 
with the increase of the thickness. 



B. /?7 sources 



After having analyzed the case of pure (3 sources, we studied (3j sources. To perform the calculations we have 
considered the 186 Re radionuclide which has the decay scheme shown in Fig. 5 01- Following Ref. (l9|. in order 
to simplify the simulation, we considered only the most relevant transitions, whose properties are given in Table llVl 
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Qec = 581.6 



122.33 . 1.69% 



7S' ; 



Qfi- = 1069.5 



2.6-10~ s % 3+ 



0.0625% 2^ 



21.54% 2+ 



70.99% 0+ 



186/-) 
76 US 



FIG. 5: Decay scheme of the 186 Re radionuclide. In actual simulations, only the /3 decays populating the ground and the first 
excited state of the 186 Os, the 7 decays shown with solid lines as well as the X-rays and IC and Auger electrons mentioned in 
the text have been considered. 



In this table, we show for X-rays, internal conversion (IC) and Auger electrons (these last two have been grouped) 
the average values of the energies and the intensities. The two fi decays considered are first forbidden (not unique). 
However, the sampling procedure [Tf3 | is the same as for allowed transitions. 



Type 


Energy [keV] 


Absolute intensity [%] 


P~ 


1069.50 


70.99 


P~ 


932.34 


21.54 


7 


767.50 


0.03 


7 


630.34 


0.03 


7 


f37.16 


9.42 


7 


122.30 


0.60 


W K X-rays 


~65.00 


6.00 


Os K X-rays 


~65.00 


3.50 


IC+Auger 


~14.00 


22.00 



TABLE IV: Characteristic radiations considered in this paper for the decay of the Re. The energies given for the f3 are 
the endpoint energies. Data have been obtained from Refs. [rH.ll9]|. 

As in the previous case, we assumed a source of R = 1.75 cm with 5 = 1 mm and an activity concentration of 
1 MBq cm -3 introduced in the inner volume of the craniopharyngioma. 

We began by comparing to the theoretical predictions based on the Loevinger and Berger approaches, our MC 
simulations, done by assuming a unique medium (water). 

For the /3 decay with _E" lax = 1.069 MeV, the mean energy is Ep — 0.362 MeV and the parameters entering in the 

Loevinger formula are c = 1.5 and v = 17.78 cm 2 g _1 . For that with E^ ax = 0.932 MeV, we used Ep = 0.309 MeV, 
c = 1.5 and v = 21.61 cm 2 g _1 . The results corresponding to -D™ w (a;) are shown in Fig. 6. The calculations have 
been done by using the absolute intensities given in Table ITVl As we can see, the simulation (dots) and the analytical 
calculation (solid lines) show a good agreement, as it happened in the case of 32 P. 

The results obtained for the remaining radiations in Table ITVl are shown in Fig. 7. Also in these results, the absolute 
intensities of each radiation are included. For photons (these are the four 7 emissions and the X-rays), the results 
of the simulations (dots) are compared with those provided by the analytical approach of Berger (solid lines). The 
parameters entering in Eq. IjAlljl . and given in Tabled have been obtained by performing cubic interpolations of 
the values quoted in Ref. [9(. It is remarkable the excellent agreement between the analytical approach and the MC 
simulation. It is also worth to notice that the 7 emissions with E 1 — 137.16 keV and the X-ray of W and Os are the 
most important contributions here, while the other emissions are more than one order of magnitude smaller. 

In Fig. 7 we also show the results the IC and Auger electrons (diamonds). These phenomena cannot be described 
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FIG. 6: Same as Fig. 2 but for the two main [3 decays of 186 Re. Only the Loevinger approach results are shown. The absolute 
intensities of both emissions are included. 
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FIG. 7: Same as Fig. 2 but for the main 7 and X-ray radiations of 186 Re. Solid curves show the result obtained from the 
Berger approach. The absolute intensities of the different emissions are included. 



by the analytical approaches considered here but, in any case, their contribution seems to be negligible. 

In Fig. 8 we show the separate contributions of the electrons, say (3 emissions, IC and Auger electrons (squares), 
and of the photons, 7 emissions and X-rays (triangles). The solid curve corresponds to the analytical calculation 
performed by summing the contributions obtained with the Loevinger approach for the two (3 emissions and with 
the Berger approach for the 7 emissions and the X-rays. As we can see it provides a rather good description of the 
simulated results. Besides, it is worth to point out that the dose delivered to the craniopharyngioma is mainly due 
to the (3 radiation, the photon contribution being two orders of magnitude smaller. On the contrary, the dose due to 
electrons becomes more than one order of magnitude smaller than the one delivered by 7 emissions and X-rays for 
x > 2 cm, that is at distances above 2 mm far from the outside surface of the craniopharyngioma. 

Finally, we have tested the effect of the different materials and interfaces in the dosimetry. First, we have evaluated 
the dose rate by assuming gel in the inner volume of the craniopharyngioma. The results corresponding to the relative 
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TABLE V: Parameters to be used in the analytical approach based on the Berger formula discussed in this paper. The medium 
considered is water. The different radiations correspond to the 7 and X-ray emissions taken into account for 186 Re. The value 
of the parameter bo is fixed to 1 for all energies. 







A v's' — vs 


(cyst) [%] 


Type 


Energy [keV] gw-ww 


gb-ww 


P~ 


1069.50 


-20.37 


-52.21 


P~ 


932.34 


-20.90 


-54.34 


7 


767.50 


-0.33 


-3.26 


7 


630.34 


-1.86 


-4.30 


7 


137.16 


-0.97 


31.55 


7 


122.30 


-1.51 


47.30 


X-rays 


~ 65.0 


-0.24 


217.38 


IC+Auger 


~ 14.0 


-20.25 


-57.89 


Full spectrum 


-20.26 


-51.19 



TABLE VI: Relative differences in %, as given by Eq. (IIOL for the different radiations of 186 Re. Here R — 1.75 cm and 
8 — 1 mm. Last row gives the results for the full spectrum of 186 Re considered in this paper. The energies given for the f3~ are 
the endpoint energies. 

differences, as given by Eq. IjlOjl . are those in the third column of Table IvTl In this table we also show the values for 
each of the radiations considered in the spectrum of the 186 Re and, in the last row, for the full spectrum (by including 
the proper absolute intensities). As we can see, the presence of the gel produces a reduction of the dose delivered to 
the cyst wall. This reduction is of ~ 20% for the f3 emissions and for the IC and Auger electrons. The reduction, no 
greater than 2%, is very small for 7 and X-ray radiations. The full effect of the gel is 20.3% reduction, as expected 
due to the dominance of the (3 radiations previously discussed. 

In a second step we have studied the effect of the material forming the cyst wall. We have checked that, as it 
happened for the 32 P radionuclide, the presence of soft tissue, instead of water, produces negligible modifications. In 
the last column of Table IvTl we show the results obtained when compact bone is assumed to constitute the cyst wall. 
As we can see, the reduction of the dose delivered to the cyst reaches more than 50% for the /?, IC and Auger electron 
emissions. However, the behavior for 7 radiations is different. While for the more energetic photons the dose reduces 
(less than 5%), it increases considerably for the less energetic ones, with an enormous enhancement (more than 200%) 
for X-rays. In any case, the dominance of the (3 emissions on the full spectrum is again remarkable. 

The average energies of the j3 emissions in this radionuclide are smaller than these of 32 P. This implies that 
Bremsstrahlung effect is less important in this case. On the other hand, we have check that the Bremsstrahlung effect 
for the 7 emissions is negligible, the dose rates obtained without Bremsstrahlung overlapping with those found in the 
full calculation. 

To finish, we compare the integral dose rates (see Eq. (|ll|l ) delivered by the two radionuclides here considered inside 
the cyst wall. In the case of a unique medium (water) and assuming the same radionuclide concentration injected in 
the inner volume of the cyst, the dose rate deposited in the wall by the 198 Re is 33% smaller than the one due to 32 P. 



12 



o 



1 + 



10 2 



10' 



10" 



Iff' 



10' : 




186 Re 

(uniform distribution; 1 MBq cm 3 ) 
R = 1.75 cm 
5 — 0.1 cm 

/T+IC+Auger 
7+X-ray 



1.5 2.0 2.5 



3.0 3.5 

x fcml 



4.0 4.5 



5.0 



FIG. 8: Same as Fig. 2 but for the full spectrum emission of 186 Re. The separate contribution of the f3 emission plus IC and 
Auger electrons (squares) and that of 7 plus X-rays (triangles) are shown. Solid curve shows the analytical result obtained 
from the Loevinger and Berger approaches. 



When we consider the cyst filled with gel and the cyst shell made of bone, the dose rate of 198 Re is 25% smaller than 
that of 32 P. 



IV. CONCLUSIONS 



In this paper we have investigated the dosimetry for radiocolloid therapy of cystic craniopharyngiomas. The 
craniopharyngioma has been described as a spherical shell with given internal radius and thickness. Different materials 
have been considered to form the colloid introduced in the craniopharyngioma and the shell representing the cyst 
wall. Explicit analytical expressions for the dose rate due to (3 and 7 emitters have been obtained on the base of the 
well known Loevinger and Berger formulae for point sources. Results for 32 P and 186 Re are quoted. 

The results obtained with these analytical expressions, valid under the assumption of a unique medium, have been 
compared with calculations performed with the Monte Carlo simulation code PENELOPE. These calculations show 
a nice agreement for both f3 and 7 emissions in the two radionuclides studied. 

By using the Monte Carlo simulation, the role of the different materials and of the interfaces present in the problem, 
has been investigated. The main conclusions we draw are the following ones. 

1. The effect of the material forming the colloid (gel in our case) produces a reduction of ~ 20% of the dose 
delivered to the cyst wall by (3 and IC and Auger electron emissions. This reduction is smaller than 2% for 7 
and X-ray emissions. 

2. There are no noticeable differences between the results obtained when soft tissue or water are considered to 
form the cyst wall. 

3. The consideration of compact bone as the material forming the cyst wall gives rise to relevant effects. For f3 
and IC and Auger electron emissions, the dose delivered to the cyst wall reduces by a factor larger than 25%. 
This effect together with that produced by the colloid makes this dose to be reduced by more than 40%. For 7 
and X-ray emissions, the presence of the bone in the cyst wall appears to be very important for energies below 
200 kcV. In this situation, the dose delivered to the cyst wall increases by more than 200% in the case of the 
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X-rays of 65 keV we have considered. For the 186 Re the dominant emissions are /3 radiations, and they obscure 
this enormous effect. 

4. All the effects quoted above, do not change with the size (inner radius) of the craniopharyngioma, but increase 
with its thickness, which appear to be the relevant parameter to consider in the dosimetry. 



APPENDIX A: FORMULAE FOR SPHERICAL UNIFORM SOURCES 

In this appendix we discuss the procedure to calculate analytically the integral in Eq. (JTJ for the point dose rates 
corresponding to the Loevinger and Berger approaches. 

In the first case, and taking into account Eq. (J2J, the dose rate per activity unit at a distance x from the center of 
the distribution can be written as 

jDp(x) = h{x) + h(x), (Al) 

where 

I t (x) = f drTi(|x-r|). (A2) 




FIG. 9: Geometrical scheme to calculate the integrals of the Loevinger and Berger formulae for an extended spherical source 
uniformly distributed. 



The key point to evaluate the corresponding integrals is to fix the center of the coordinate system on the target 
point. The scheme is represented in Fig. 9. In the case of the integral Ii(x), one has to consider that T\{x) has a 
range t (see Eq. |J3J|). Thus, the only contribution to the integral comes from points inside the source distribution 
and sited at a distance s < t from the target point. Therefore, only the shadowed volume shown in Fig. 9 contributes 
to Ii (x) and we have 

P pt fOrtx&x /•27T 

h{x) = / dsTi(s) = / dss 2 ! 1 !^) / dtfsinfl / d(j> 

JV Jx-R JO Jo 

ft /-COS #max 

= 2tt / dss 2 Ti(s) / d(-cosfl) 

Jx-R Jo 

= 2tt f dss 2 (l-cos0 max ) Ti(s). (A3) 

Jx-R 

The angle 8 max is the value reached by 9 when |x — s| = R and then 

x 2 + s 2 — R 2 

cos 6 max = . (A4) 

2xs 

Considering the above equations, the integral (|A3|) becomes 

7T r* 

h{x) = Q(x- R)G(R + t-x) dss [{x - s) 2 - R 2 ] Ti(s). (A5) 

x Jx-R, 
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Here, the Q(x — R) factor indicates that the calculation is valid for points outside the distribution, while the 8(i?+< — x) 
factor points out the fact that this contribution cancels for x values larger than R + t. The result of the integration is 



h(x) = 



nBt 
~2x~ 



12x- llt + 4(t-R) exp 
2{R 2 -x 2 )\n' ' 



t + R- x 



3i? 2 + 2Rx - 5x 2 



x-R 



e(x-R)e(R + t-x) 



Note that, in the limit x = R this contribution is 



h(R) 



nBt 
2R 



{AR [3 - exp(l)] - t [11 - 4exp(l)]} 



(A6) 



(A7) 



To calculate the contribution l2(x), a similar procedure is followed. Now, the full source distribution contributes 
and this implies that in Eq. 1|A5|I we must change the upper limit in the integral to x + R. Thus one has 



h{x) 



7T f X+R 

- dss [(x- s) 2 - R 2 ] T 2 (s) 

x Jx-R 



2nBi 



■ exp 



z — R — x 



2R 

(R + z) + {R- z) exp ( — 



By substituting Eqs. l|A6j) and (|A8|> into Eq. IjAlj) we have finally 



nBt 
~2x 



2-kBz 2 



12x- lit + A(t - R) exp 



t + R-x 
t 



(A8) 



(A9) 



3R 2 + 2Rx - bx 2 + 2{R 2 - x 2 ) In 



x-R 



■ exp 



z — R — x 



2R 

(R + z) + (R- z) exp ( — 



Q{R + t-x) 
@{x - R) . 



For 7 emitters, the integral in Eq. JJ| is evaluated for the point dose rate given by the Berger approach ©. In this 
case, J-y(x) has no range and the calculation follows the same steps as those for the T 2 term of the Loevinger formula. 
Thus one can write 



1 f it r x +R 

-rDj(x) = / drJ 7 (|x-r|) = / dss [(x - s) 2 - R 2 ] J 7 (s). 

Substituting Eqs. © and @ here one has 



-D^x) 



1 N+3 

1 j=0 



(A10) 



(All) 



where 



rj( x ) = 



min(i,2) 

^ ai(x)bj-i [i^ 1 , 

i=max(0 ,j-N—l) 



(A12) 



with 



Oi(x) = 




(A13) 



and 



K s (x) 



x+R 



ds 1 exp(— /is) . 



(A14) 



x-R 
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If j = the integral Kq{x) is given by ^(| 



K Q {x) = Ex iji{x - R)) - E x (fi(x + R)) (A15) 

= ln£± l + g< = iTi f : x _ 

n—l 



where E\ (x) is the exponential integral ^(j ■ For j ^ we have 01 



A,(x) = -L [ 7 (i, fi{x + R))--y (j, nix - R))} , j?0, (A16) 



with 



7 (n, z) = (n - 1)! P(n, z) = (n - 1)! i 1 - exp(-z) T\ \ ( Al7 ) 
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